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Abstract 

A long-standing open question asks for the minimum number of vectors needed to form 
an unextendible product basis in a given bipartite or multipartite Hilbert space. A partial 
solution was found by Alon and Lovasz in 2001, but since then only a few other cases have 
been solved. We solve all remaining bipartite cases, as well as a large family of multipartite 
cases. 
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1. Introduction 

We use "ket" notation to represent unit vectors in complex Euclidean space, which 
represent pure quantum states. We say that a state \v) G C dl g) • • • <g) C dp is a product state 
if we can write it in the form 

\v) = \vi) <g> • • • ® | Up) with \vj) G <£ d] V j. 

If no such decomposition of \v) exists, we say that it is entangled. 

An unextendible product basis (UPB), introduced in [BDM + 99, DMS + 03], is a set of 
mutually orthogonal product states such that every state in the orthogonal complement of 
their span is entangled. That is, a set S C C dl (g) • • • <g) C dp is a UPB if and only if: 

(a) \v) is a product state for all \v) G S; 

(b) (v\w) = for all \v) ^ \w) G S; and 

(c) for all product states \z) S, there exists \v) G S such that (v\z) ^ 0. 



Email addresses: ciLenkensliiii@gmail.com (Jianxin Chen), nathaniel.johnston@uwaterloo.ca 

(Nathaniel Johnston) 



Many applications of UPBs have been found, including the construction of bound en- 
tangled states [BDM + 99, LSM11, Skoll] and indecomposible positive maps [TerOl]. Un- 
extendible product bases have also been shown to give rise to Bell inequalities without a 
quantum violation [ASH + 11] and feature the interesting property of nonlocality without 
entanglement [BDF + 99] - that is, they can not be perfectly distinguished by local quantum 
operations and classical communication, even though they contain no entanglement. 

One of the oldest questions concerning UPBs asks for their minimum size - the smallest 
number of vectors that form a UPB. A trivial lower bound of /jv(di, . . . , d p ) := ~ 
1) + 1 was immediately noted in [BDM + 99], and it was also shown that this lower bound 
holds in many cases. In fact, if we remove the orthogonality condition (b) in the definition 
of a UPB, then this lower bound is always attained [Bha06]. However, when condition (b) 
is present (which is the case we consider), the problem seems to be more difficult. It was 
shown in [AL01] that there are cases when the lower bound /jv(^i, • • • , d p ) is not attained, 
and furthermore it was determined exactly for which values of d±, . . . , d p this is the case. 

The minimal size of UPBs has only been determined in a handful of cases when the 
trivial lower bound is not attained. Our main contribution to this problem is a solution for 
a large family of systems, which includes all bipartite (i.e., p — 2) systems as a special case. 
We review known partial answers to this question in Section 2 before presenting our main 
results. We also briefly introduce our proof technique, which is based on factorizations of 
the complete graph and tools from algebraic geometry. In Section 3 we present the full proof 
of our main result, and in Section 4 we prove a related result that answers the question for 
some additional tripartite (i.e., p = 3) systems. As our proofs are largely non-constructive, 
we discuss how to explicitly construct UPBs of the minimal size in Section 5. In order 
to illustrate these procedures, we provide MATLAB code and several examples in small 
dimensions. We close in Section 6 with a brief discussion of related questions that remain 
open. 

2. Minimum Size of Unextendible Product Bases 

A long-standing open question asks: given the dimensions d%, rf 2 , . . . , d p , what is the 
minimum possible number of vectors in a UPB S C <C dl <E>- ■ -cg)C dp ? We define f m (di, . . . , d p ) 
to be this minimum value. Here we briefly list all partial answers to this question that are 
known to us: 

(1) If min(<i 1 , d 2 ) = 2 then f m (di, d 2 ) = did 2 . 

(2) f m (dx, . . . ,d p ) equals the trivial lower bound /at(g?i, • • • , d p ) := Y?j=i{dj — 1) + 1 if and 
only if (1) doesn't hold and either /jv(rfi, . . . ,d p ) is even or all d/s are odd (or both). 

(3) lip = 2 (mod 4) then / m (2, 2, . . . , 2) = p + 2 [Fcn06]. 

(4) / m (4, 4) = 8 [Pcd02], f m (2, 2, 3) = 6, f m (2, 2, 5) = 8, f m (2, 2, 2, 2) = 6, f m (2, 2, 2, 4) = 8, 
and / m (2, 2, 2, 2, 5) = 10 [Fcn06]. 
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Notice in particular that in all known cases except for (1) and (2) above, we have 
f m (di, . . . , dp) = fw(di, . . . , dp) + 1. Our main result shows that this is a fairly general phe- 
nomenon, and there is a rather large class of multipartite systems for which f m {d\, . . . , d p ) = 
f N (di, ...,d p ) + l. 

Theorem 1. Let 2 < d\ < d 2 < • • • < d p be integers for which neither (1) nor (2) above 
hold, and d p — 1 > Y^Z\{dj — 1) > 3. Then f m (d\, . . . , d p ) = /jv(c?i, . . . , d p ) + 1. 

This result generalizes the known results / m (4, 4) = f m {2, 2, 2, 4) = 8 and / m (2, 2, 2, 2, 5) = 
10. More importantly, Theorem 1 solves all remaining bipartite cases. The following corol- 
lary follows immediately from Theorem 1 and the known results (1) and (2). 

Corollary 2. In the p = 2 case, we have the following characterization: 

(a) If mm(di,d 2 ) = 2 then f m (di,d 2 ) = d^d 2 . 

(b) If di, d 2 > 4 are even then f m (di, d 2 ) = d\ + d 2 . 

(c) If neither (a) nor (b) hold then f m (di,d 2 ) = d\ + d 2 — 1 . 

Note that the only way that the condition Y^j=i(^j ~~ -0 — 3 of Theorem 1 can fail is 
if either p = 2 and di < 3 (in which case the size of the minimum UPB is given by either 
condition (1) or (2) above), or if p = 3 and di = d 2 = 2. Thus, the theorem says nothing 
about the values of f m (2,2,2k + 1) for k > 3. The following result solves half of these 
remaining cases: 

Theorem 3. Let k be a positive integer. Then f m (2, 2, 4k + 1) = /at(2, 2, 4k + 1) + 1. 

The next two sections are devoted to proving Theorems 1 and 3. Our techniques are 
very much inspired by the orthogonality graph and graph factorization methods introduced 
and developed in [BDM+99, DMS+03, Fen06]. 

We illustrate the basic ideas briefly here, before providing the full proofs in the next 
sections. In all cases, condition (2) above tells us that f m (dx,...,d p ) > /jv(^i, • • • , d p ) + 1, 
so it suffices to construct a UPB of size /W(<ii, . . . , d p ) + 1. For now, assume for simplicity 
that we are trying to construct a UPB S := {\v ), . . . , |t> 5 ), \w ), . . . , \w 5 )} 6 C 6 ® C 6 . For 
all < j < 5, we can write \vj) = \v^) <g) \vj ) and \wj) = \u)j) <S> {w^}, and we know that 
any two of these product states are orthogonal on at least one of their subsystems (e.g., for 
all i 7^ j we have either (v\ ) = or (v\ \vj ) = 0, or both). 

In order to easily keep track of these orthogonality conditions, we use orthogonality 
graphs. All graphs that we consider will be simple graphs. The orthogonality graph of a set 
of vectors {|0o), • • • , is the graph (V, E) defined by letting V := {0o, • • • , 4>k-i\ be a 

set of k vertices, and E := {((pi, 4>j) : (4>i\(fij) = 0} be a set of edges so that two vertices are 
adjacent if and only if the corresponding vectors are orthogonal. 

Our goal now is to construct sets of vectors := {\v^}, . . . , \vjp), \ w^}, . . . , Iw^)} G 
C 6 and := {\v^}, . . . , \v^), I^q 2 '')' • • • j l^s 2 " 1 )} e ^ 6 so ^ na ^ their orthogonality graphs 
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(V, Ei) and (V, ^2) satisfy (V, Ei U -E2) = K12, the complete graph on 12 vertices, since this 
ensures that any two distinct elements of S are orthogonal. For convenience, in the future if 
we have two graphs G\ := (V, E\) and G2 '■— (V, E 2 ), then we define G\ UG2 '■= (V, E\ UE^)- 
We begin by considering the graph D 6 := (V,Ei), which is constructed by taking two 
copies of Kq and pairwise joining vertices, as in Figure 1. We also consider a second graph, 
C^o := (V, E 2 ), which we define simply to be the complement of D 6 (i.e., an edge is in E 2 
if and only if it is not in Ei), as in Figure 2. 



Wq 




V 3 



Figure 1: The graph Z?6,o used in the construction of UPBs in which d p — 1 = Y^=i(dj~ 1) = 5. 



V Vi v 2 v 3 v 4 v 5 




Wq Wi W2 W3 W4 W5 



Figure 2: The graph C^o, which is the complement of D§ o- 

We now must show that we can construct two sets of 12 vectors in C 6 that have or- 
thogonality graphs D e o and Cq^. For C^o, this certainly seems like it should be possible, 
since each vector lives in C 6 and is orthogonal to exactly 5 other vectors - we use tools from 
algebraic geometry to make this intuition rigorous. 

For Dqq, the idea is to define \vj ) := \j) to be computational basis states for all 
< j < 6. Then if we want the orthogonality graph of to equal Dq q, we need to find 
states {|wj 2 ' ) )} that form an orthonormal basis of C 6 and satisfy (w^lj) = for all j. In 
other words, we need to find a unitary matrix that has zeroes on its diagonal. We again 
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use methods from algebraic geometry to prove that many such unitary matrices exist, and 
furthermore that generic sets of vectors with these orthogonality graphs lead to product 
bases that are indeed unextendible. The ideas presented above generalize without much 
difficulty to all cases considered by Theorem 1. 

3. Proof of Theorem 1 

We now present the main body of the proof of Theorem 1. We make use of the two 
technical lemmas 4 and 5 throughout the proof, whose statements and proofs we leave to 
the end of the section. 

Define the quantity 6 := ((d p — 1) — 5^=i( £ £j — l))/2 > 0, which we can assume is an 
integer since condition (2) of Section 2 solves the case when it is a half-integer. Define Dd Vyb to 
be the graph on2(d p -b) = f N (d l7 d p ) + l vertices V := {v , . . . , v dp „ b -i, w , . . . , w dp _ 6 _i} 
such that Vi is adjacent to Vj and Wi is adjacent to Wj for all i ^ j, and Vi is adjacent to Wj 
if and only if j — i = s (mod (d p — b)) for some < s < b. In the d p = 6, b = case, this 
graph is depicted in Figure 1. The graph Dj t i looks the same, except it has 6 additional 
edges: (v , w x ), (v t , w 2 ),..., (v 5 , w ). 

Our first goal is to construct a set of 2(d p — b) vectors 

s {p) ■= {\v^h . . . , 1^}, K w >, • • • , G 

such that the orthogonality graph of is D dp b . To this end, let \v^) := \j) be standard 
basis states for all < j < d p — b. We then need to find states {|w^)}^L b 1 that form 

an orthonormal set and have the additional property that ((j + £) (mod d)^^) = for all 
< £ < b. We choose these states to be the (normalized) columns of a matrix U described 
by Lemma 4, which is proved at the end of this section. Importantly, note that any d p + 1 
distinct elements of span all of C dp . To see this, suppose we choose r + 1 |uk^)'s and 

dp — r |fj p ^)'s from (for some b < r < d p — b) . It then suffices to show that every r x (r + 1) 
submatrix of U has full rank r. When 6 + 2 <r < d p — b — 2, this claim follows immediately 
from condition (Hi) of Lemma 4. Similarly, condition (in) says that every (b + 2) x (b + 2) 
submatrix of U has full rank, so every (b+1) x (b + 2) submatrix of U must have rank 6+1, 
which proves the r = 6 + 1 case. The r = d p — b — 1 case is similar. The last remaining case 
is r = 6, in which case it is not true that every r x (r + 1) submatrix of U has full rank r. 
Nonetheless, because all d p — 6 \vj )'s have been chosen, it suffices to just show that every 
6x6 submatrix contained entirely within the bottom 6 rows of U is nonsingular, which is 
guaranteed by (iv) of Lemma 4. This completes the proof that any d p + 1 distinct elements 
of span all of CV 

Now we consider the complement of the graph Dd p ,b, which we denote Cd p ,b- This is a 
graph on the same set of 2(d p — b) vertices V, but (v i: Vj) and (u>j, Wj) are never edges in C^b, 
and (vi, Wj) is an edge in Cd v ,b if and only if j — i = £ (mod (d p — 6)) for some 6 < £ < d p — b. 
For example, in the d p = 6, 6 = case, this graph is depicted in Figure 2. The graph CV,i 
looks the same, but without the edges (vo, u>i), (v±, W2), ■ ■ ■ , (V5, wo). 
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We now split C dp) b into p—1 spanning subgraphs - one for each of the first p—1 systems. 
Let Cd p ,b,r, s (r,s > 1 and s + r < d p — b) be the subgraph of Cd p ,b that contains the edge 
(vi, Wj) if and only if j — i = t (mod (d p — b)) for some s < £ < s + r. For example, the 
graphs (^6,0,2,1 an d ^6,0,3,3 are depicted in Figure 3. To help keep in mind which indices 
represent what, we can think of r as the regularity of the graph (i.e., Cd p ,b,r,s is r-regular 
- every vertex has degree r), and s as a shift (since, for example, the smallest index j 
such that (vo,Wj) is an edge in Cd P) b,r,s is j = s). It is straightforward to see that we have 

Cd p ,b, ri , Sl U • ■ - UCd p ,b,r k ,s k = C dp ,b whenever Ylj=i r j = dp -26-1 and si = 1, Sj = sj-i+rj-i 
for 1 < j < k. We choose k = p — 1 and rj = dj — 1 for all 1 < j < p — 1. 



W 1>1 ^2 ^3 W 4 V 5 D Di D 2 D 3 f 4 f 5 




IOq U7i W 2 W 3 Wi W 5 W W\ W 2 W 3 W4 W 5 



Figure 3: The graph C^o, decomposed into the union of C6,o,2,i and C6,o,3,3- 

Our goal now, for each 1 < j < p, is to find a set of 2(d p — b) vectors S^) whose 
orthogonality graph is Cd p ,b,d —i, a • We choose S^' to be the (normalized) columns of the 
matrix V described by Lemma 5 (with q = d p — b, r = dj — 1, and s = Sj). Not only do 
these vectors have the desired orthogonality graph, but property (ii) of the lemma ensures 
that any dj distinct elements of (1 < j < p) span all of C dj . 

Since we have C dpA d 1 -i, Sl U • • • U C dpAdp „ 1 -i, Sp _ 1 U D dp , b = C dp>b U D dp>b = K 2 ( dp -b), the 
complete graph on 2(d p — b) vertices, it follows that the set of product states 




are mutually orthonormal. To see that this set is unextendible, recall that any d p + 1 of the 
vectors in span all of C dp and any dj of the vectors in span all of C dj for 1 < j < p. 
Thus, any product state not in S can be orthogonal to at most d p elements of S on the p-th 
subsystem and at most dj — 1 elements of S on the j-th subsystem for 1 < j < p, and thus 
it can be orthogonal to a total of at most YljZi(dj ~ ^) + d p = f N (di, . . . , d p ) elements of 
S. However, there are f N (dx, . . . ,d p ) + 1 vectors in S, so no product state is orthogonal to 
them all, which completes the proof. 

We now present and prove the two technical lemmas that we made use of earlier in this 
section. 

Lemma 4. For all integers b > and d > 4 + 2b, there exists a matrix U G Md t d-b such 
that the following four conditions hold: 

(i) the columns of U form an orthonormal set; 
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(ii) ((j + t) (mod (d - b))\U\j) = for all < j < d - b and < £ < b; 

(Hi) allrxr sub-matrices ofU are nonsingular whenever r G {6+2,6+3, . . . , d— b}\{d— b— 1}; 
and 

(iv) allbxb submatrices ofU that are contained entirely within the bottom b rows ofU are 
nonsingular. 

Proof. For any U G M d ^ d _ b) we write 

/ x 1A + iy 1A x h2 + ij/i.2 • • • x^d-t, + ij/i,d_6 \ 
_ | ^2,1 + ^2,i ^2,2 + ^2,2 • • • a?2,«i-6 + ij/ 2 ,d-& 

\ £d,i + ij/d,i x dj2 + iy d;2 • • • + iy d ,d-6 / 

where y p>q G K. for all 1 < p < d, 1 < q < d — b. 

Then conditions (i) and (ii) can be rewritten as follows: 

(i 1 ) x ij+e)imod ( d -b)),j = y(j+£)(mod (d-b)),j = for all 1 < j < d - b and < i < b, 
d 

(ii') ( x k P ~ iyk P ){xkq + iykq) = 0Wl<p<q<d-b, 
k=i 

We are mainly interested in the real variety Z C R 2cf ( d_fe ) characterized by conditions 
{%') and (ii'). We also look into the real variety ■^j 1 j 2 ...vjij 2 -> — M? d ( d ~ b ^ characterized 
by polynomials of x p ^s and 2/ p , ? 's arising from the constraint det(C/j 1 j 2 ...j T . J1 j 2 ...j T .) = 0, 
where U ili2 ... ir j 1 j 2 ...j r is the r x r submatrix of [/ formed by rows ii , i 2 , ■ ■ ■ ,i r and columns 

Jl,j2, • • -,jr- 

Our aim is to show that 

^\(( U Ziii2-ir,jij2-jr) U ( U % i\ii-ib,hh-3b) ) 7^ 

re{i)+2,6+3,...,d-6}\{d-b-l} d-i>+l<»i <•••<»;,<<! 

l<ij < ■■<i r <d l<j 1 <---<j b <d-b 
l<jl<---<jr<d-b 

Note that if ■^i 1 i 2 ...i T .,j 1 j 2 ...> fl Z is a proper subset of Z, then -^i 1 i 2 ...i r ,j 1 j 2 ...> H Z has zero 
measure in Z. In order to show that ^iii 2 ...i r ,jij 2 -> H Z is a proper subset of Z, we first need 
some simple definitions. A permutation matrix is a square binary matrix that has exactly one 
entry 1 in each row and each column and 0's elsewhere. Similarly, a column-permutation 
matrix is a (not necessarily square) binary matrix that has exactly one entry 1 in each 
column, at most one entry 1 in each row, and 0's elsewhere. Clearly, every permutation 
matrix W G M d>d has full rank d and every column-permutation matrix W G M^d~b has 
rank d — b. 

Let W G M djd _ b be the matrix with a 1 in the ((j + b + l)(mod (d — b)), j)-entry for all 
1 < j < d — b, and 0's elsewhere. The top (d — b)x(d — b) submatrix of W is a permutation 
matrix, so it is nonsingular. Therefore the corresponding point in M. 2d ( d ~ b ^ that represents 
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W lies in Z\ Z 12 ...(d-b),i2...(d-b), which implies that Zi2...(d-&),i2...(d-6) H Z is a proper subset 
of Z. Hence Zi2...(d-6),i2...(d-&) H Z has measure zero in Z. 

Similarly, for any 1 < i\ < ■ ■ ■ < id-b < d, we can always move some appropriately chosen 
nonzero rows of the above W to its lower 6 x (d — b) submatrix such that the corresponding 
point in R 2d( - d - b ) lies in some Z\Z hi2 _ i(d _ b))12 ...(d-b)- Hence, ZnZ ili2 _ i(d _ b))12 ...(d-b) has measure 
zero in Z. 

Next, for any fixed 1 < i\ < • ■ ■ < i r < d and 1 < ji < ■ • • < j r < d — b, where 
6 + 2 < r < d — b — 2, if there is some column-permutation matrix W G Md^d-b satisfying 
condition (ii) such that W / iij 2 ...ir,jij2...jr * s a permutation matrix, then the point in R 2d ( d_fe ) 
that represents W lies in Z \ (Zi2...(d-&),i2...(<z-&) U ^»it2...irjii2...>) > ^ follows that if such a 
matrix W exists then Z D (-Zi2...(d-6),i2...(d-&) U ^hi2-i r ,hh-3r) a ^ so has measure zero in Z. 

We can construct G Md^d-b to satisfy these requirements as follows. We choose 
the submatrix formed by rows i\,% 2 , . . . ,i r and columns ji,j 2 , . . . ,j r to be a permutation 
submatrix and the submatrix formed by rows {1,2, ... ,d} \ {ii,i 2 , . . . ,i r } and columns 
{1, 2, . . . , d — b} \ {ji,j 2 , , ■ ■ ■ ,j r } to be a column-permutation submatrix. All other entries 
are chosen to be zero. 

To see that the above two submatrices of W can be chosen as desired, note that condition 
(ii) forces (b + l)(d — b) entries of W G Md,d-b to be zero. For any given r x r submatrix 
■ir,hj2-jr formed by rows i\,i%, . . . ,i r and columns ji,j 2 , . . . ,j r , where b + 2 < r < 
d — b — 2, there are at most (6+1) preset zero entries in each column j t (1 < t < r). 
So we can first fix a nonpreset entry to be 1 and then seek an (r — 1) x (r — 1) 

permutation submatrix Wi a ,„i T) j 2 ...j r . This observation implies that we only need to prove 
the existence of (6 + 2) x (6 + 2) permutation submatrix formed by any (6 + 2) rows and 
(6 + 2) columns that satisfy condition (ii). We also require the existence of (6 + 2) x 2 
column-permutation submatrix formed by any (6 + 2) rows and any 2 columns that satisfy 
condition (ii), but this claim is obvious. Hence, the only thing left is to prove is that, given 
fixed 1 < i\ < ■ • • < ib+ 2 < d and 1 < j\ < • • • < jb+ 2 < d — 6, there exists W G Md t d-b 
satisfying condition (ii) such that Vt^ 1 i 2 „.j 6+2 j 1 j 2 ...j I>+2 is a permutation matrix. 

To see that such a matrix W exists, note that there is no x x y submatrix of W with 
x + y = b + 3 that is forced entirely to equal zero by condition (ii). Thus the same is true 
of all (6 + 2) x (6 + 2) submatrices of W. Thus there is no (6 + 2) x (6 + 2) submatrix of 
W that satisfies condition (ii) of [HS93, Theorem 3.1] with m = n = 6 + 2, r = 6 + l. By 
then using condition (Hi) of [HS93, Theorem 3.1] and Lemma 3.7 of the same paper, we see 
that every (6 + 2) x (6 + 2) submatrix of W can be made nonsingular by a suitable choice of 
the entries that are unaffected by condition (ii). If we write the determinant as the sum of 
products of entries of the (6 + 2) x (6 + 2) submatrix, then at least one of the terms in the 
sum can be made nonzero - the permutation matrix we want is the one that corresponds to 
this term in the sum. 

We have shown that Z fl Zi 1 i 2 ...i r ,jij2-jT h &s measure zero in Z for any r G {6 + 2, 6 + 
3, . . . , d — b} \ {d — b — 1}, 1 < it < ■ ■ • < i r < d and 1 < ji < ■ • • < j r < d — b. 

For any d — b-\-l<ii <■•■<% <d and 1 < ji < ■ ■ ■ < jt < d — 6, we choose suitable 
W G Md,d-b suc h that its 6x6 submatrix formed by rows i\, i 2 , . . . , % and columns ji,j 2 , ■ ■ ■ ,jb 
is a permutation matrix and all other entries are fixed to be zero. Its corresponding point 
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in R 2d ( d fe ) lies in Z \ Zi 1 i 2 ...i b ,j 1 j 2 ...j b - Hence, Z n ^i^.-.i^hn-h ^ as measure zero in Z. 
This implies that 

^fl (( U Zi 1 i 2 ...i r j 1 j 2 ...j r )u( Zi\ii-i h ,nn...jS) 

re{b+2,b+3,...,d-b}\{d-b-l} d-b+l<i 1 <---<i b <d 
l<i!<---<i r <d l<j 1 <...<i i) <d-6 
l<Jl<'"<Jr<d-6 

also has measure zero in Z since the union of a finite number of measure zero subsets is 
again a measure zero subset. In particular, it follows that 

^\(( U Ziii2-ir,jij2-jr) U ( [J ^iii2-i b ,hh-3b) ) 7^ 0- 

re{i)+2,6+3,...,d-6}\{d-b-l} d-6+l<ij <---<i b <d 

l<il<-<i r <d l<j 1 <---<j b <d-b 
l<j 1 <---<j r <d-b 

This implies the existence of matrix U G M^d-b with properties (i), (ii), (in) and (if). 
Indeed, a generic matrix [/ G M^a-b satisfying (i) and (ii) will also satisfy conditions (iii) 
and (if). □ 

Lemma 5. Let q,r,s be positive integers satisfying q > r + s. Then there exists a matrix 
V G M r+ i j2 q such that: 

(i) for all < j < q and s < £ < s + r , V\j) is orthogonal to V\(j + £)(mod q) + q); and 

(ii) every (r + 1) x (r + 1) submatrix ofV is nonsingular. 

Proof. For any V G M r+li2 ,, we write 



V = 



( X lt i + X li2 + ZJ/1,2 • • • Xi t2q + iyi,2q \ 

X2,l + m,l X 2 ,2 + iy2,2 ■ ■ ■ X 2 ,2q + ^2,2<? 

\ Xr+1,1 + Wr+1,1 X r+ i t 2 + ij/r+1,2 ' ' ' 2V+l,2g + Wr+l,2q J 



where Xij,yij G R for all 1 < % < r + 1, 1 < j < 2q. 
Condition (i) implies that 

U\V*V\(j + e)(modq) + q) = (1) 

for all < j < q and s < £ < s + r. Equivalently, 

r+l 

y]( x k,j+l - Wk,j+l)(Xk,(j+£)(mod q)+q+l + Wk,(j+e)(mod q)+q+l) = ( 2 ) 
k=l 

for all < j < q and s < £ < s + r. 

So every matrix V G M r+lj2(J satisfying condition (i) corresponds to a solution of the 
above polynomial system. Let Z C R 4 ( r+1 )9 be the real variety that is characterized by this 
polynomial system. We also look into the real variety Zj^ 2 ...j r+1 C ]R 4 ( r+1 )9 characterized 

9 



by polynomials of the Xijs and yijs arising from the condition det{Vj 1 j 2 ,„j r+1 ) = 0, where 
Vj 1 j 2 ,,,j r+1 is the (r + 1) x (r + 1) matrix formed by columns ji,j 2 , ■ ■ ■ ,jr+i of V. 
To show that 

l<h<-<jr+l<2q 

it suffices to prove that ZnZ w - 2 „j r+1 is a proper subset of Z for any 1 < j 1 < ■ ■ ■ < j r+1 < 2q. 

To see why this is sufficient, note that if this is the case then Z D Zj 1 j 2 ...j r+1 has zero measure 

in Z for all 1 < j% < • • • < j r +i < 2q. A finite union of measure zero subsets is again a 

measure zero subset, which implies Z D ( |J ^hh-Or+i) nas measure zero in Z. 

i<ii<-<i,+i<2 g 

The lemma then follows straightforwardly, since a generic matrix in M r+ i y2q that satisfies 
condition (i) will also satisfy condition (ii). 

To complete the proof, we thus just need to show that, for any fixed 1 < j\ < ■ ■ ■ < 
jr+i < 2g, there exists some V G M r+1)2(J such that, for all < j < q and s < £ < s + r, 
V'\j) is orthogonal to V'\(j + £)(mod q) + q) and the (r + 1) x (r + 1) submatrix formed by 
columns ji,j2, ■ ■ ■ ,jr+i of V is nonsingular. 

To this end, we first choose the ji,j2, ■ ■ ■ , jr+i-th column vectors of V to be the set of 
computational basis states in (r + l)-dimensional space. Thus the submatrix formed by these 
columns is nonsingular. We then fill the other columns with suitable computational basis 
states one by one. To fill the k-th column vector, note that it is orthogonal to r other column 
vectors by condition (i). If r' of these r column vectors have already been determined, then 
we can choose the k-th column vector to be any computational vector that is orthogonal to 
the r' < r already-determined column vectors. By repeating this procedure, we eventually 
fill the entire matrix V so that it satisfies condition (i) and satisfies condition (ii) for the 
particular submatrix formed by columns jx, j 2 , ■ ■ ■ ,j r +i- D 

4. Proof of Theorem 3 

Our proof of Theorem 3 is similar in style to that of Theorem 1 . As it is already known 
that / m (2, 2, Ak + 1) > f N (2, 2, Ak + 1) + 1, it suffices to find a UPB of size f N (2, 2, 4k + 1) + 1. 
To this end, we begin by presenting a graph that leads to a product basis of the desired size, 
and then we use algebraic geometry techniques to show that many of these product bases 
are unextendible. 

We begin by defining the graph Y^+i on /at(2, 2,4/c + 1) + 1 = 4k + 4 vertices V : = 
{v o, • • • , v 2 k+i, Wo, ■ ■ ■ , w 2 k+i} such that Vi (wi) is adjacent to Vj (wj) if and only if j — i G 
{1,-1} (mod {2k + 1)), and v , is adjacent to Wj if and only if % = j (this graph is sometimes 
called the prism graph on 4k + 4 vertices). In the k = 2 case, this graph is depicted in 
Figure 4. 

We now construct a set of vectors in C 2 ®C 2 whose orthogonality graph is Y^+i- To this 
end, we follow the notation of [Fen06] and let {\bj), l&j - )} 2 ^ 1 De distinct orthonormal bases 
of C 2 (i.e., (bj\bf) = for all j, but 1(6^)1, |(^|^)| $ {0, 1} whenever i ^ j). We 
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^5 ^0 




^3 V 2 



Figure 4: The graph Y\2 used in the construction of a minimal UPB in C 2 ® C 2 ® C 9 . 
then define the vectors |fj \ G C 2 <8> C 2 as follows: 

\ V 2j' 2) ) \ b j) ® l«2i+l) := ® l^(2j+2)(mod(2fc+2))) 

l W 2i 2) ) : = ® l™2j+l) : = ® l&(2i+3)(mod(2fc+2))), 

for all < j < k. It is straightforward to verify that the set of vectors 

sm ^ {1^) i«fiS>, k m ) I«SS» 

has orthogonality graph Y^ + ^. Furthermore, any product vector \z) G C 2 <E> C 2 can be 
orthogonal to at most 2 elements of S^ 1 ' 2 ) on the first subsystem, and at most 1 element of 
gP-, 2 ) on the second subsystem. Thus any nonzero product vector can be orthogonal to at 
most 3 elements of S^ 1,2 \ 

We now consider the complement of the graph Y ik+4i , which we denote X 4fc+4 . This is a 
graph on the same set of 4fc + 4 vertices V, but (v^Vj) and (wi,Wj) are edges in X 4fc+4 if 
and only if j — i £ {1, —1} (mod (2k + 1)), and (vi,Wj) is an edge in X 4 fc +4 if and only if 
i ^ j. Our goal is to show that there exists a set 

of 4k + 4 vectors with orthogonality graph X 4fc+4 . To this end, we let be the columns 
of the matrix W described by Lemma 6. It is clear from conditions (i), (ii) and (Hi) of the 
lemma that these vectors have the desired orthogonality graph. Furthermore, condition (iv) 
guarantees that any nonzero state \z) G C 4fe+1 can be orthogonal to at most 4k vectors in 

Because X 4 fc +4 Ul 4 fc +4 = K^+i, the complete graph on 4k + 4 vertices, the set of vectors 

s ■= {|«?' 2) > ® k 3) >, • • • , |t&S> ® k 1,2) > ® k 3) >, ■ ■ ■ , l«4!£) ® 

is a product basis of C 2 £§> C 2 ® C 4fc+1 . To see that S is unextendible, simply recall that 
any product state can be orthogonal to at most 3 elements of S on the first two subsystems 
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and at most 4k elements of S on the third subsystem, for a total of at most 4k + 3 states 
total. Thus no nonzero product state can be orthogonal to all 4k + 4 elements of S, so 
unextendibility follows and the proof is complete. 

We finish this section with the statement and proof of Lemma 6, which played a key role 
in the above proof of Theorem 3. 

Lemma 6. Let k>l. There exist matrices W^,W^ € M^+i^k+i such that if we define 
the block matrix W := [W( v ), W( w )] £ Mt& . | |/, . . ; then: 

(i) W( v )\i) is orthogonal to W( v )\j) whenever j —i^ {0, 1,2k + 1} (mod (2k + 2)); 

(ii) W(w)\i) is orthogonal to W(w)\j) whenever j — i {0 ,1,2k + 1} (mod (2k + 2)); 

(Hi) W( v )\i) is orthogonal to W( w )\j) whenever j ^ i; and 

(iv) every (4k + 1) x (4k + 1) submatrix of W is nonsingular. 

Proof. Define d := 4k + 1 for simplicity. Similar to the proof of Lemma 5, to show the 
existence of W , it suffices to prove the following claim: 

For any 1 < j 1 < ■ ■ ■ < j d < d + 3, there exists some W G M djd+3 that satisfies the 
orthogonality conditions (i), (ii) and (in), and the d x d submatrix formed by columns 
ji, 32, ■ ■ ■ , jd is nonsingular. 

We first choose the ji, J2, • • • , J<rth column vectors to be the set of computational basis 
states in d- dimensional space. The submatrix formed by these columns is clearly nonsingular. 
We then fill other column vectors with suitable computational basis states one by one. To 
fill the k-th column vector, note that conditions (i), (ii) and (Hi) force it to be orthogonal to 
(d — 1) other columns of W. If we assume that r of these (d— 1) columns have already been 
specified, then we can choose the k-th column vector to be any computational basis state 
that is orthogonal to the r < d — 1 already-specified columns. By repeating this procedure, 
we eventually fill the entire matrix W so that it satisfies conditions (i), (ii), and (Hi), and 
satisfies condition (iv) for the particular submatrix formed by columns . . . ,j d - □ 

5. Explicit Construction of UPBs 

As the proofs of Lemmas 4, 5, and 6 are non-constructive, it is perhaps not immediately 
clear how to produce explicit UPBs of the size indicated by Theorems 1 and 3, even though 
we know they exist. We now address this problem and demonstrate how to construct explicit 
UPBs of the desired size in small dimensions. Code that implements the techniques described 
in this section, and thus constructs minimal UPBs, can be downloaded from [Johl2]. 

For Lemma 5, we recall that a generic matrix satisfying condition (i) of the lemma will 
also satisfy condition (ii). Thus, one way to construct matrices satisfying both requirements 
is to randomly generate its first q columns, then generate its last q columns according to 
the orthogonality requirement (i), and finally check to make sure that the resulting matrix 
satisfies condition (ii). All of these steps are straightforward to implement numerically in 
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software such as MATLAB. For example, the following matrix ^5,3,2 satisfies the conditions 
of the lemma in the <? = 5, r = 3, s = 2 case: 



5,3,2 



3 3 3 11-1-1 2-2 

2 1 1 2 3 -2 -2 -1 

2 1111 5-3-4 2 1 

22322 2 1 3 



The procedure for Lemma 6 is similar - just randomly generate columns that satisfy 
conditions (i), (ii) and (in), and then with probability 1 the resulting matrix will also 
satisfy condition (iv). More specifically, randomly generate the first two columns of the 
matrix, then generate the remaining columns according to the orthogonality conditions, and 
finally check to make sure that condition (iv) is satisfied. As before, this procedure is simple 
to perform numerically. The following matrix W2 is an explicit example in the k = 2 (i.e., 
d = 9) case: 



Wo 



2 


1 


1 


5 


1 


2 








1 





6 


6 


2 


-1 


-1 


-3 


1 


2 








-1 





6 


6 


1 


2 





-4 


4 


1 








2 


1 


3 


3 


2 


2 








-4 


-1 








2 


1 


27 


6 


2 


2 














-1 


-1 


-4 








-4 


1 


2 














1 


2 


-4 








-2 


2 


2 




















3 


-1 


-20 


-14 


1 


2 























-1 


21 


23 


2 


2 


























-31 


-12 



Unfortunately, constructing matrices that satisfy the constraints of Lemma 4 seems to 
be a bit more difficult in practice. We can begin by specifying (d — h)(d — b — l)/2 of the 
non-zero entries and then using the orthogonality condition (i) to solve for the remaining 
non-zero entries. In general, however, this leads to a system of (d — b)(d — b — l)/2 linear 
and quadratic equations in (d — b)(d — b — l)/2 variables, which quickly becomes infeasible 
to solve as d grows. Furthermore, the solution will generally be significantly messier than 
we saw when dealing with Lemmas 5 and 6. Nonetheless, such a system of equations can be 
solved (even by hand) when d is small enough. For example, in the d = 6, b = 1 case, if we 
set 



U e ,i 






u 5 
1 

u s 

1 



ill 




u 7 

1 
1 



1 

«3 




Ug 
1 



U-2 
1 

u 6 


1 





U4 
1 
1 




and then solve for the variables {%}}=i, we find that the unique real solution is obtained 
when Uio is the unique real root (approximately equal to 1.6445) of the polynomial 3wf — 
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2u\ Q — 3u 
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3, and the remaining u/s are given by 

w 4 = u 10 (u w - 2)/ (2uio - 3) m 7 = -uio 

u 5 = -l-u 10 u 8 = u w -l 

u 6 = 1/(1 + u w ) it 9 = 1/(1 - « 10 ) 



mi = (uio - 2)/(l - uio) 
u 2 = {uw - l)/(uio - 2) 
it 3 = (3 - 2«io)/(uio - 2 



It is a simple calculation in MATLAB to verify that U^i also satisfies conditions (Hi) and (iv) 
of Lemma 4, and thus leads to an unextendible product basis. 

In the 6 = case, however, we can simplify things a bit further by analytically construct- 
ing a matrix U^q that satisfies conditions (i) and (ii) of the lemma. We start by building 
the eigendecomposition of Ud,o and then we argue that it must satisfy the condition (ii). 

Let the first two eigenvalues of Ud,o be Ai := 1 and A2 := —1, with respective eigenvectors 



1 T 



and |f 2 ) :- 



1 T 



V^2' V23=2' ' ' ' ' s/2d^2 



It is easily-verified that \v 1) and \v%) are orthogonal, and a simple calculation reveals that 



+ A 2 |f2)(w 2 | 





1 



/d^T 
1 

/d^T 



/d^T 



/d-1 





/d-1 




/d-1 




Now let := e 2m / k be a primitive fe-th root of unity. Let the remaining d — 2 eigenvalues 



of Ud,Q be the (d — 2)-th roots of unity: A J+2 
corresponding eigenvectors as follows: 



d-2 



for j 



1,2, 



d — 2. Define the 



' J j+2 



fd-i 



2j 
^d-l' 



j(d-2) 
d-1 



In other words, [^3) , ^4), . . . , \vd)\ is the (d — 1) x (d — 1) Fourier matrix, with its leftmost 
column removed and a row of zeroes added to the top. The fact that {|fj)}j=i forms an 
orthonormal basis of C d comes from the fact that the (d— 1) x (d— 1) Fourier matrix is unitary. 
It follows that £7^0 := X^/=i ^j\ v j)( v j\ * s unitary as well. Furthermore, a straightforward 
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(albeit tedious) calculation reveals that if k, £ > 1 then 



3=1 
d-2 

_ 1 \^ , ,j(d-l+(fc-Q(d-2)) 
~ d - 1 2^^-1)^-2) 

(d-2)(d-l+(fc-*)(<i-2)) _ 1 
U7 (d-l)(d-2) l 

~ 7~, 777 i-«i-(fc-fl(d-2)N' 
(a - W - w (d _ 1)(d _ 2) J 

where we summed the geometric series in the final line above. This quantity equals zero 
if and only if (d - 2)(d - 1 + (A; - i)(d - 2)) is a multiple of (d - l){d - 2) (i.e., if and 
only if (k — £)(d — 2) is a multiple of d — 1). However, because < \k — £\ < d — 1 and 
— 2 is coprime to d — 1, this happens if and only if k = £. This verifies that U^o indeed 
has zeroes down its diagonal, as desired, and also shows that it is a candidate to satisfy 
condition (m) of the lemma (since in the 6 = case, any matrix satisfying condition (Hi) 
must have all of its off-diagonal entries be non-zero). We have numerically verified that U^o 
satisfies condition (Hi) (and it trivially satisfies condition (iv)) for 4 < d < 19 in MATLAB. 

By using all of these techniques, together with the techniques used in the proof of The- 
orem 1, we can explicitly construct UPBs in all of the dimensions we have discussed. For 
example, a minimal UPB of size 10 in C 4 Cg> C 6 is given by the set 

® l^ 2) >, • • • , \vP) ® \v?), k 1} > ® ki 2) ), . . k 1} > ® |4 2) )}» 

where: 

• the \vj )'s and |w^)'s are (in order) the normalized columns of ^5,3,2 above; 

• \ V W) = \j) e C 6 for < j < 4; and 

• |wj 2 ' ) ) is the (j + l)-th column of C/ 6j i above. 

6. Outlook 

We have shown that, in many cases, the minimum size of a UPB does not exceed the 
trivial lower bound by more than 1. In fact, there is currently no known case in which 
f m (di, . . . , d p ) > /jv(c?i> • • • , d p ) + 1. It could be the case that this never happens, or it could 
be the case that we aren't aware of any such cases yet because it is very difficult to prove 
non-trivial lower bounds on f m (di, . . . , d p ). 

Some particularly interesting cases of the minimal UPB question that remain open are: 
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(1) d 1 = d 2 = 2, d 3 = 4Jfe-l: It was shown in [Fen06] that f m (2, 2, 4Ac— 1) = /jv(2, 2, 4fc— 1)+1 
when fe = 1, but the proof technique does not seem to generalize straightforwardly to 
the k > 2 case. 

(2) p = 4k and d\ = ■ ■ ■ = d p = 2: All other cases with d\ = ■ ■ ■ = d p = 2 have been solved - 
it is known that / m (2, . . . , 2) = f N (2, . . . , 2) if pis odd and / m (2, . . . , 2) = /y(2, . . . , 2)+l 
if p = 2 (mod 4). Furthermore, it is known that / m (2, 2, 2, 2) = /W(2, 2, 2, 2) + 1, but 
again the proof technique does not obviously generalize to the k > 2 case. 

(3) d\ = 3, g?2 = d% = 4: Excluding the open case (1) above, this is now the smallest 
unsolved tripartite case. 

Finally, it was noted in [AL01] that whenever f m (di, . . . , d p ) = /at(^i, • • • , d p ), a minimal 
UPB can be constructed using only real vectors. The same is true in all cases in which 
we have proved / m (c?i, . . . , d p ) — /at(^i, • • • , d p ) + 1. Restricting to the real field cuts the 
number of variables in half and also cuts the number of polynomial constraints in half, but 
has no significant effect on any of our proofs. 

Acknowledgements. Thanks are extended to John Watrous for a helpful discussion. J.C. 
was supported by NSERC, UTS-AMSS Joint Research Laboratory for Quantum Computa- 
tion and Quantum Information Processing and NSF of China (Grant No. 61179030). N.J. 
was supported by an NSERC Postdoctoral Fellowship. 



References 



[AL01 
[ASH+11 
[BDF+99 

[BDM+99 
[Bha06 
[DMS+03 

[Fen06 
[HS93 
[Johl2 
[LSM11 
[Ped02 



N. Alon and L. Lovasz. Unextendible product bases. J. Combinatorial Theory, Ser. A, 95:169- 
179, 2001. 

R. Augusiak, J. Stasinska, C. Hadley, J. K. Korbicz, M. Lewenstein, and A. Acm. Bell inequalities 
with no quantum violation and unextendible product bases. Phys. Rev. Lett., 107:070401, 2011. 
C. H. Bennett, D. P. DiVincenzo, C. A. Fuchs, T. Mor, E. Rains, P. W. Shor, J. A. Smolin, 
and W. K. Wootters. Quantum nonlocality without entanglement. Phys. Rev. A, 59:1070-1091, 
1999. 

C. H. Bennett, D. P. DiVincenzo, T. Mor, P. W. Shor, J. A. Smolin, and B. M. Terhal. Unex- 
tendible product bases and bound entanglement. Phys. Rev. Lett., 82:5385-5388, 1999. 

B. Bhat. A completely entangled subspace of maximal dimension. Int. J. Quantum Inf., 4:325- 
330, 2006. 

D. P. DiVincenzo, T. Mor, P. W. Shor, J. A. Smolin, and B. M. Terhal. Unextendible product 
bases, uncompletable product bases and bound entanglement. Commun. Math. Phys., 238:379- 
410, 2003. 

K. Feng. Unextendible product bases and 1-factorization of complete graphs. Discrete Appl. 
Math., 154:942-949, 2006. 

D. Herschkowitz and H. Schneider. Ranks of zero patterns and sign patterns. Linear Multilinear 
Algebra, 34:3-19, 1993. 

N. Johnston. Code for computing some unextendible product bases. Published electronically at 
http : //www . nj olmston . ca/publications/minimum-upbs/ code/, 2012. 

J. M. Leinaas, P. 0. Sollid, and J. Myrheim. Unextendible product bases and extremal density 
matrices with positive partial transpose. E-print: arXiv:1104.1318 [quant-ph], 2011. 
T.B. Pedersen. Characteristics of unextendible product bases. Master's thesis, Aarhus Univer- 
sitet, Datalogisk Institut, 2002. 

16 



[Skoll] L. Skowronek. Three-by-three bound entanglement with general unextendible product bases. J. 

Math. Phys., 52:122202, 2011. 
[TerOl] B. M. Terhal. A family of indecomposable positive linear maps based on entangled quantum 

states. Linear Algebra AppL, 323:61-73, 2001. 



17 



